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Abstract 

A model of the universe as a very large white hole provides a useful alternative 
inhomogeneous theory to pit against the homogeneous standard FLRW big bang 
models. The white hole would have to be sufficiently large that we can fit comfortably 
inside the event horizon at the present time, so that the inhomogeneities of space-time 
are not in contradiction with current observational limits. A specific Lemaitre-Tolman 
model of a spherically symmetric non-rotating white hole with a few adjustable 
parameters is investigated. Comparison of calculated anisotropy in the Hubble flow 
and the CMB against observational limits constrain the parameter space. A Copernican 
principle would require that we are not too near the centre of the white hole. As an 
additional constraint this predicts a value of Do between 0.9999 and 1. 
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Introduction 

If recent observations are taken at face value and a zero cosmological constant is 
assumed, then the universe is open with 0.2 < .Q < 0.7 [1]. The standard Friedmann- 
Lemaitre-Robertson- Walker (FLRW) cosmology gives a consistent homogeneous 
model for such an open universe, but in the instant after the big bang it describes an 
infinite, very nearly flat, expanding universe. The horizon problem is especially bad in 
this case no matter how much inflation there may be. The cosmological principle that 
the universe is homogeneous on scales larger than the observable universe is therefore 
more difficult to accept for open universes than it was for closed. Although Q ^ 1 is 
certainly not yet excluded, this may be a good time to re-examine the simplest open 
inhomogeneous cosmological models which do not suffer from the infinite horizon 
problem. 

The isotropy of the Cosmic Microwave Background (CMB) is taken as very good 
evidence for homogeneity on the scale of the observable universe from about 100 
million lyrs to 10 billion lyrs, but we have no evidence for homogeneity far beyond the 
observable horizon. One possible inhomogeneous model of the universe would be a 
spherically symmetric, non-rotating white hole. There is a common misconception that 
this would be inconsistent with observation because a white hole has tidal forces which 
are not observed in the cosmos. This is not the case because tidal forces can be absent 
within the expanding dust cloud of a white hole. 

An example which demonstrates this is the model of star which collapses to form a 
black hole in the solution of Oppenheimer and Snyder [2] . The star is a sphere of dust 
of uniform density with zero pressure. Within the sphere the space-time is FLRW and 
outside it is Schwarzschild. The universe can be modelled as a white hole which is the 
time reversal of a very large version of this solution to the equations of general 
relativity. From the point of view of an observer whose observable universe lies within 
the sphere there is no way to distinguish between this white hole model and a standard 
FLRW cosmology. 

Inhomogeneous models 

If the Oppenheimer-Snyder solution was our final cosmological model it would be 
rather hopeless. There would be no possibility of testable predictions distinct from 
those of FLRW without waiting an unknown time for the edge of the dust cloud to 
come into view. However, a more physical solution would be expected to have a non- 
uniform matter density with no sharp edge. In such a cosmology deviations from 
homogeneity might be observable but could be kept within empirical limits. 

Spherically symmetric, non-rotating, dust filled cosmologies were first studied by 
Lemaitre who found a general solution in 1933 [3]. They were further developed by 
Tolman in 1934 [4]. The same Lemaitre-Tolman (L-T) model was reconsidered by 
Datt in 1938 [5], Bondi in 1947 [6], Omer in 1949 [7] and many others after (see 
Krasinski's bibliographical review [8]). The purpose of this paper is to consider 
whether a particular example of a L-T white hole is a viable model of the universe. 
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If the universe fitted a standard FLRW dust filled cosmology with zero cosmological 
constant, astronomers would have only two parameters to measure before they could 
work out the entire shape of the universe (e.g. Cl and Ho corresponding to density and 
expansion rate). The L-T model has an infinite number of parameters because the dust 
mass density and radial velocity can be arbitrary functions of radius. However, we 
could select suitable functions characterised by two parameters characteristic of size, 
or density and expansion rate. Astronomers would then have to measure these two 
plus the present age of the universe and our distance and direction (two angles) from 
the centre of spherical symmetry, a total of six parameters. Since it is already difficult 
to measure Do and Ho it is not likely to be easy to determine these extra parameters 
which can only appear as anisotropy and inhomogeneity. 

A direct approach is to plot the radial and angular dependence of the Hubble 
expansion. Suitable data from Type la supernovae is now being accumulated and it 
should soon be possible to set tight limits on the Hubble anisotropy. Another 
possibility would be to study the CMB radiation. The situation is less discouraging if 
we can believe that the dipole anisotropy has a measurable contribution from very 
large scale inhomogeneity rather than just from the influence a local great attractor. 
The residual dipole could be determined by subtracting the velocity of the Local Group 
from our apparent velocity relative to the CMB. This would give the direction to the 
centre of the white hole model and one other parameter. Higher order moments in the 
CMB might provide further input if they are not due to statistical anisotropy on the 
surface of last scattering as conventional wisdom supposes. Thus there is some hope 
that the parameters of the white hole model can be measured if suitable anisotropics 
are not masked by local effects. 

The objective for this paper is therefore to try to construct a non-rotating white hole 
model of the universe consistent with observations and in which the CMB dipole and 
Hubble anisotropy is explained with or without a great attractor. Relations between 
anisotropics will be predicted. 

There has been previous work on the effects of inhomogeneity on the CMB using the 
L-T models, notably by Raine and Thomas in 1981 [9], Paczynski and Piran in 1990 
[10], Kurki-Suonio and Liang in 1992 [11], Panek in 1992 [12], and Arnau, Fyullana 
and Saez in 1994 [13]. Some of these have considered the effect of an inhomogeneity 
situated within the observable universe so that the surface of last scattering is not 
perturbed. Here I will consider the effects of inhomogeneity which extends through the 
observable horizon. 
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The LemaTtre-Tolman Model 



The Lemaitre-Tolman Model (sometimes called the Tolman-Bondi model) is the 
general, spherically symmetric, pressureless dust solution for the gravitational field. We 
shall consider only the case of zero cosmological constant. In co-ordinates co-moving 
with the dust, the metric takes the form, 

ds 2 = dt 2 - [1 + 2£(r)r 1 R, 2 (t, r)dr 2 - R 2 (t, r)[d0 2 + sin 2 (6)d(p 2 ] 



E(r) is an arbitrary function which can be loosely interpreted as total energy per unit 
mass of the dust shell at r. The area and circumference of the mass shell at r are that of 
a sphere of radius R(r,t) but the actual distance to the centre of symmetry may be more 
or less than that, depending on whether E(r) is positive or negative . The radial 
parameter r can be rescaled using any monotonic function without changing the form 
of the metric 

R(t,r) must satisfy the equation, 

R, 2 = 2E + 2M(r)/ R 

where M( r) is another arbitrary function which can be interpreted as K times the total 
dust mass within the sphere out to r. The form of this equation is identical to the 
energy equation for a test particle travelling radially out from a sphere of mass M(r)//C. 
It will escape if E>0 or collapse back if E<0. The matter-density is given by, 

P(r,t) = 



K R R 



'r 



(K = 87TG) 

For an exploding white hole out of which all matter will eventually escape, E(r) is 
positive everywhere and the solution is given by, 



^ Af 



V J 



cych(;c) is the hyperbolic cycloid function defined by 

cych(sinh T]-T]) = cosh 77 - 1 
cych'(^) = s]l + 2/ cych(x) 

cych(.*;) = (f)M + * f -lio* 2 + tf(* f ) 
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to(r) is an arbitrary function describing the location of the initial singularity, and is 
known as the bang time. At large time this solution describes a dust sphere with shells 
expanding at velocity, 



Eliminating The Bang-Time 

In order to further simplify the model we eliminate the arbitrary bang-time function 
to(r) by setting it to zero. This step can be justified on physical grounds. A non- 
constant bang-time function would mean that the singularity was not simultaneous in 
co-ordinates co-moving with the flow of matter. Provided the singularity is space-like 
there will always be some reference frame in which it is simultaneous. This is a 
reference frame with cosmo logical time defined as the longest proper- time on any 
timelike curve which extends back from each event to the singularity. So a statement 
equivalent to setting zero bang-time is that the initial flow of matter is stationary in a 
synchronous reference frame with cosmological time. In a small region near the 
singularity the universe is very close to being homogeneous. If we assume that the flow 
of matter is determined in a causal fashion from the physics of the early universe, then 
the large scale deviations from homogeneity cannot have any significant effect. By the 
symmetry of the situation in a small causal domain it follows that the flow cannot be 
biased in any direction, i.e. we expect a zero bang-time. 

A second argument for a zero bang-time is an application of a variational principle near 
a singularity. The field equations of general relativity are derived from the principle of 
least action applied over an internal region of space-time. A stationary point of the 
action is found when the fields are varied over the region but fixed on a boundary. In 
the presence of singularities this is an incomplete answer. It should also be necessary to 
consider what happens over a region which includes a singularity along the boundary. 

Consider again the case of a free dust filled space-time and vary the dust flow field 



The dust field is varied by displacing each particle of dust by an amount given by an 
arbitrary vector field \f . The corresponding change in the flow field is 





Sp»=(p v bv- p vb v ). v 
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So varying the action gives 

= 2\(p- l p fl y y (p'b"-p l 'b')4^gd > x 

-2\p- l p fl (pV-p> , b*)dS„ 

s 

The first term determines the well known equation of motion for the dust. The second 
is a boundary term at the singularity which is zero for all \f iff 

p"p v dS v -p"p^dS x =0 

dS A is a vector orthogonal to the singularity. The equation requires that the dust flow is 
parallel to this vector. In other words, that the singularity is simultaneous in co-moving 
co-ordinates. 

There is one other significance of a zero bang-time which is worth mentioning. The 
Weyl curvature is smaller near the singularity. The Weyl curvature in the L-T model 
has components given by, 

C = — K 

R 3 6 

Near the singularity the L-T solution is approximately, 



R(t,r) = -M(ry(t-t (r)) 



, , 2 M, r 6 M 
p(r,t) = -^- = - 

K R R, K R 



1 + 2 M f o, 



C = -2 



M 



R 3 



M t , r 



So to first order the Weyl curvature vanishes when the bang-time is a constant function 
of r. In that case it still diverges as t tends to zero to, but less fast than the Ricci 
Curvature. It only vanishes completely in the homogeneous FLRW limiting case. 

Shape of the White Hole 



To go much further we must specify a shape for the white hole model of the big bang 
by specifying functions for E(r) and M(r). Recall that the metric is invariant under 
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reparameterisations of r so there is really only one arbitrary function to specify. The 
homogeneous open FLRW case is, 

E(r) = \E r 2 
M(r) = M r 3 

Because of the flexibility to rescale r, this is really a one parameter model. The 
Oppenheimer-Snyder model is a two parameter solution which is the same as above for 
r<l but with M(r) constant for r>l. 

For a finite mass white hole in general, M(r) must be a monotonic function which is 
bounded as r tends to infinity. If E(r) has an r-squared behaviour for small r then M(r) 
must be r-cubed. If the mass distribution is to be smooth then these can be smooth 
functions with E(r) = E(-r), M(r) = -M(-r). This still leaves much freedom of choice 
but to be specific we take, 

£(r) = ^£ sinh 2 (r) 
M(r) = M tanh 3 r 

This can be written in other forms by reparameterisations, e.g., 

E(r) = {r 2 

M r 3 
M(r)= 



,3 

L ^2 



(E + ry 

The first form is the more convenient and will be used here. 

Matter Distribution 



R(r,t) = (iM(r)t 2 f + 
P(r,t) = 



81 



E(r)f 



500 JI//V \ 

V M ( r ) J 



- 27 E(rf t 2 +. 



M(r) 



K R 2 R, r 



4 (32 ^ 3ME, r +M, r E 

\375J 1 4 



3m KM'M^t 



3 



The first term of R(r,t) remains finite as r tends to infinity and the first term of the 
density is independent of r. This means that near the singularity almost all the matter of 
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the white hole is concentrated within a finite shrinking volume of almost uniform 
density. It is a good approximation to the Oppenheimer-Synder solution for a 
collapsing dust sphere. As time progresses the mass spreads out and the edge becomes 
less distinct. 

This alleviates the horizon problems which would apply to the open FLRW cosmology 
which has an infinite size near the singularity. It also means that a universe based on 
this model will be approximately homogeneous until near the edge of the dust sphere. 
An observer whose past light cone passed outside the dust sphere would see a large 
hole in the CMB. Since the CMB appears very isotropic to us we can assume that we 
are well within the sphere. 



Hubble Parameters 

In a non-homogeneous universe the Hubble expansion can have different rates in 
different directions. Observation suggests that the Hubble parameter is nearly 
independent of direction although variations of up to 10% might be tolerated. 
Measuring it accurately is difficult because it is small and probably masked by local 
inhomogeneity. Setting limits on its size could nevertheless be a very useful constraint 
on our possible position in the white hole. 

For the white hole cosmo logical model we define a radial Hubble parameter H R and a 
transverse Hubble parameter H T . We define them from the time dependence of the 
metric components. The Hubble anisotropy H A is the relative difference of the two. 



H , — 



H A = 



__ R, rt _ 2 ( 4 V 3ME,-EM, r 




M M. t 



(3ME,-2EM, r )t 
M*M,„ 



The Hubble anisotropy is small for small t or small r so the low observed anisotropy 
indicates that we are either near the beginning of an especially large white hole, or near 
the central axis. Using the conjectured functions for the mass and energy, Ha can be 
plotted against r and s = (E 1/2 /M 1/s ) t 1/s , See figures la and lb. (These plots use exact 
calculations rather than the first order formulae given above.) 

For small s, 



H A = 



20 (2f 



sinh 2 (2r)/+. 
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Omega 

Observations also indicate the ratio of the local average matter density to the critical 
density required to cause a recollapse. 

H = l(H R + H T ) 

Kp 



Q = 



3H 2 

' (6ME, r +EM, r )t 



= 1- 



V250) 



'r 



This can also be plotted against r and s, see figures 2a and 2b. 
For small s, 



l-£l = r cosh 2 (r)[2cosh(2r) + 3]s 2 +. . . 

1()(2) ! 

If observations confirm that Q is not close to one then s must be quite large. For 
example, if Q is 0.5 then s is about 1.0 near the central axis of the white hole. A small 
Hubble anisotropy means that we are near the central axis. For example, a Hubble 
anisotropy of 5% corresponds to r of about 0.6. 

If known together, the Hubble anisotropy and Q determine r and s. The Hubble 

3 2 

parameters also determine t so the ratio E '/M is then also known. 



Cosmic Microwave Background Radiation 

Further data can be taken from measurements of the CMB. It is well known that a 
CMB dipole is observed and that this is interpreted as the proper motion of our solar 
system. Two contributions are suspected; our motion round the galaxy and the motion 
of our galaxy towards the great attractor. The Milky Way galaxy appears to be moving 
at 600 km/sec (0.002c) towards Virgo relative to the CMB. This motion could be 
interpreted as a sign of large scale anisotropy rather than a local motion. More 
conservatively, we can only say that the residual dipole after local effects have been 
removed is not likely to be much bigger than this. Observations of the Hubble flow 
should be able to determine the motion of our local group. This can then be subtracted 
from the apparent motion relative to the CMB to yield the true CMB dipole in co- 
moving co-ordinates. 

The CMB also has measurable quadrapole and higher order components. These are 
traditionally attributed to temperature variations at the "surface of last scattering" 
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when radiation decoupled from matter in the early universe. They could also be due to 
inhomogeneity in the curvature of space-time either at the surface of last scattering or 
afterwards. Since these anisotropics are very small (about 1 part in 100,000) they 
could place strong constraints on the white hole model. 

In order to predict the dipole effect from the white hole model we must calculate the 
CMB in directions towards and away from the centre as seen on a co-moving geodesic 
(r = const) . The white hole model will predict a zero quadrapole moment because of 
the rotational symmetry about any axis which passes through the centre. To find the 
octopole and higher order anisotropics we must calculate the redshift as a function of 
the angle <J between the direction of observation and the direction to the centre. 

The spectrum of the observed radiation will also depend on the shape and temperature 
of the surface of last scattering which may in turn depend on unknown physics of the 
early universe. In the standard model it is found that after last scattering the 
temperature is related to matter density as Tocf}' 3 '. If we assume that the energy 
available is predominantly from processes which reached equilibrium then this relation 
will also hold independently of expansion rates at the time. It will be almost universal 
even in an inhomogeneous universe. 

As we have already seen, shortly after the big bang the matter density was very 
uniform. The moment of last scattering was very early and we expect the present size 
of observable universe to be relatively small compared to the extent of the white hole. 
It is therefore safe to assume that the CMB temperature was homogeneous at the 
outset. First order contributions to the anisotropy will come from the subsequent 
evolution of the geometry of the universe. 

By the spherical symmetry we know that the dipole will be orientated along the radial 
direction out from the centre of the white hole. Light-like geodesies in the radial 
direction are given by, 



dt 
dr 



= + 




R 



If St is the time difference along two subsequent geodesies in the same direction then 




with the integral being along the light-curve. 
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This can also be expressed as a partial differential equation for lp = log(<^)> 

dt dr dt 

To get the dipole it is necessary to solve this in both directions and take the difference. 
This deals with the (7=0 and (J=TT directions. The more general case for any direction 
is a little more difficult because the light-like geodesic equations must be used. The 
redshift function can be shown to satisfy this partial differential equation, 

RG + -zj— R cos o — GR, t sin 2 a - RG, t cos 2 a = 
at or 

{R, r sin <7 — sin <7 cos c(RG, t -GR, t )} 

The solution can be found in the form of a series, 

y/(r, t) = y/ log(0 + £ yr(r, a)p 

i 

i = 1,2, • • • 

Any term which does not vanish for small t has to be dependent of r and <J as a 
boundary condition. Expansions for G and R are known, 



Gir^^G^' 

i 

7?(r,0 = £^(r)r fi 



i = 1,2, 



Collecting terms of the lowest order, which is t m , gives, 



Vo^i - f G X R X sin 2 a - 1 G X R X cos 2 a = 



v = ! 



Taking higher order terms in sequence we get, 
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v,=o 



= — sin C7 + — COS C7 
R x G 1 



1 



f 



\ 



\ R 1 G \J 



1 

2 



^2 G : 



V / 3=-V2>r 



COSG 



cos(2<7) 



etc. 



The octopole is the coefficient of cos(2cr), 



M 8 = 



80(2) 5 



sinh 2 (2r) l s ,/ 



and the dipole is the coefficient of cos(cr), 



2 A /l + J E sinh / (r) 



cosh 5 (r)sinh(r)^ 3 



Contributions to the octopole moment appear earlier than contributions to the dipole 
moment. This suggests that the residual dipole moment should be small compared to 
the octopole which is 10~ 5 . The only hope of a larger dipole is if E is small or if r is 
small, but if Eo is too small the edge of the mass sphere would fall inside the observable 
universe. The observed dipole of 4x10" is most likely due to local motion rather than 
large scale inhomogeneity. 

Not surprisingly there is a direct relationship between the octopole CMB anisotropy 
and the Hubble anisotropy. 



This means that there is very little hope of observing a contribution to Hubble 
anisotropy from large scale inhomogeneity. 



The Copernican Principle 

Another consideration is worth mentioning even though it is statistical. On the 
assumption that there is no special reason that we should be near the centre of the 
white hole, it is very unlikely that we should happen to be in the small volume which 
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has low values of r. We can make the ansatz that the a priori probability distribution 
for our location is given by the normalised mass density. This will be reasonable 
provided galaxies of the right sort develop everywhere with a density proportional to 
mass density. It is a big assumption, but even if it is not true there is no reason to 
suppose that the right type of galaxy should form near the centre so the ansatz may be 
robust. 

The proportion of matter within the sphere up to r is tanh 3 r. This means that we can 
say that r > 0.4 with 95% confidence. The median value is r = 1.1. This is as good as 
most co smo logical observations at present but such statistical constraints need to be 
interpreted with great care. 

If we combine this with the observed limits on the CMB anisotropy, it is clear that we 
must live in the small s region near the singularity. It also allows us to place limits on 
the value of Q : 1 > O « 0.9999. It is then unlikely that observation within the 
observable universe will be able to distinguish between closed or open co smo logical 
models. 

The only escape from this conclusion would be an explanation for why we should live 
very close to the centre of the white hole. As an example for illustrative purposes only, 
we could consider the effects of an electric charge spread over the matter of the 
universe. This would generate electric fields everywhere except near the centre and 
could make other regions inhabitable. By an anthropic argument we could then 
understand why we find ourselves near the centre. 
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